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This paper analyzes the onset characteristics of a miniature thermoacoustic Stirling heat engine using the ther¬ 
modynamic analysis method. The governing equations of components are reduced from the basic thermodynamic 
relations and the linear thermoacoustic theory. By solving the governing equation group numerically, the oscilla¬ 
tion frequencies and onset temperatures are obtained. The dependences of the kinds of working gas, the length of 
resonator tube, the diameter of resonator tube, on the oscillation frequency are calculated. Meanwhile, the influ¬ 
ences of hydraulic radius and mean pressure on the onset temperature for different working gas are also presented. 
The calculation results indicate that there exists an optimal dimensionless hydraulic radius to obtain the lowest 
onset temperature, whose value lies in the range of 0.30-0.35 for different working gases. Furthermore, the am¬ 
plitude and phase relationship of pressures and volume flows are analyzed in the time-domain. Some experiments 
have been performed to validate the calculations. The calculation results agree well with the experimental values. 
Finally, an error analysis is made, giving the reasons that cause the errors of theoretical calculations. 
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Introduction 

Thermoacoustic engine (TE) converts heat into acous¬ 
tic power and has the advantage of no moving parts and 
thus promises a long life. The TEs have been applied in 
some field such as cryogenic refrigeration and electric 
generation. In standing-wave TE [1], the energy conver¬ 
sion between the gas and solid through irreversible ther¬ 
mal contacts and the theoretical thermal efficiency is less 
than 20%. The thermoacoustic Stirling heat engine 
(TASHE), firstly proposed by Backhaus and Swift [2], 
employs an inherently reversible thermodynamic circle 
and the thermal efficiency can reach up to 30%. The lin¬ 
ear thermoacoustic theory [3] is successful in explaining 


the conversion from thermal energy to acoustic energy. 
The DeltaEc, developed by Ward and Swift based on 
linear thermoacoustic theory, is a powerful tool for opti¬ 
mizing and designing TEs. However, some nonlinear 
phenomenon, such as onset characteristics and mode 
transition, are still not fully understood. So far, few re¬ 
search works involve the onset characteristics of TE. Lai 
[4] studied the onset temperature of a standing wave TE 
based on the thermoacoustic network theory, however, 
the method cannot be applied directly to the TASHE for 
the absence of practical geometrical closed end in the 
loop; Benavides [5] predicted the onset temperature by 
taking the output work of gas circle to be zero as the 
critical onset conditions, however, the coupling of oper- 
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ating frequency and onset temperature has not been taken 
into account in his model; Hu [6] proposed the two-port 
network model and obtained the startup criteria for TE by 
using Nyquist instability criterion. Recently, de Waele [7] 
proposed a thermodynamic model to reveal the onset 
conditions in TASHEs, focusing on the onset temperature 
and oscillation frequency. In this model, some assump¬ 
tions have been made for simplicity. The regenerator is 
treated to be ideal (zero void volume), which indicates 
the compliance effect of regenerator is neglected. Mean¬ 
while, the acoustic field in the regenerator is assumed to 
be a pure traveling wave, meaning the oscillating pres¬ 
sure and volume flow are in phase. In an actual TASHE, 
however, the acoustic field in the regenerator is neither a 
pure standing wave nor a pure traveling wave but a trav¬ 
eling-standing wave [8]. The performance of TASHE, 
including the onset characteristics, is strongly dependent 
on the acoustic field. The working conditions (such as 
working gas or mean pressure) and regenerator’s geomet¬ 
rical parameters (such as porosity or hydraulic radius) 
have great influence on the acoustic field. Therefore, a 
more practical thermodynamic model needs to be devel¬ 
oped. In this paper, a similar model is built with a more 
realistic treatment of regenerator based on de Waele’s 
research. Numerical calculations and experiments are 
carried out to investigate the onset characteristics of a 
miniature TASHE. 


Inertance 







Regenerator 


Pulse tube 


JE 


i h 


I 1 ' 1 


— Buffer 


Pb 


Resonator 





Fig-1 Schematic drawing of the miniature TASHE: (a) before 
thermodynamic analogy conversion; (b) after thermo¬ 
dynamic analogy conversion. 


Governing equations 

The miniature TASHE is schematically shown in Fig. 
1(a). The main part of the miniature TASHE consists of 
cold heat exchanger (CHX), regenerator, heating block, 
pulse tube and ambient heat exchanger (AHX); together 
with feedback tube, inertance tube and 180° bend cavity, 
they form a closed loop. The resonator connects to the 
loop with a three way tube. The dimensions of compo¬ 
nents in the loop are much smaller than the wavelength 
of oscillating working gas, so they can be modeled as 
lumped elements. The void volumes of CHX, AHX and 
heating block are small enough to be negligible. The in¬ 
ertance tube is modeled as a solid piston with a mass M t 
equal to the mass of working gas contained in the iner¬ 
tance tube; only the compliance effect is considered in 
the feedback tube, 180° bend tube, pulse tube and only 
the inertance effect is considered in the resonator. The 
resonator is modelled as a cylinder with volume V R (av¬ 
erage length L r ), closed by a resonating piston with mass 
M r . The value of Mr is equal to the mass of the gas con¬ 
tained in the resonator in front of the piston. Let the basic 
angular resonance frequency be equal to the resonator 
frequency of the mass-gas spring system [2, 8], it has 
L R =2L ac /n. In this way, the TASHE can be simplified to a 
thermodynamic model as shown in Fig. 1(b). 


Table 1 Main dimensions of the miniature TASHE 


System parameter 

Sym 

Value 

Diameter of resonator 

d r 

0.016 m 

Length of resonator 

L a c 

0.35 m 

Volume of cavity 

v b 

5.4><10“ 4 m 

Diameter of regenerator 

D r 

0.021 m 

Length of regenerator 

L r 

0.02 m 

Length of pulse tube 

L t 

0.06 m 

Diameter of pulse tube 

D t 

0.022 m 

Volume of space d 

v d0 

3.84xl0 _5 m 

Length of inertance tube 

Li 

0.1245 m 

Diameter of inertance tub 

D, 

0.022 m 

Volume of space c 

V c0 

8.5xl0“ 5 m : 

Ambient temperature 

T c 

300 K 


The usual complex notation for the time-oscillating 
pressure in the system is: 


p=p„> + s p= p,„ + Re 




where p m is the mean pressure, dp is the time-dependent 
oscillating pressure, which is given by the real part of the 
product of complex amplitude p and e uot . Other vari¬ 
ables (such as volume flow rate U and volume V) have 
the same form as the oscillating pressure p. It is assumed 
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that the time-dependent oscillating variables are much 
smaller than the mean variables. In the following analysis, 
we will focus on the time-dependent oscillating variables 
and all the dynamic equations are in time-domain. 

The equation of volume d is given by: 

v = ^K l+ Ko_^P l (2) 

d dt yp 0 dt 

The dynamic equation of mass of the inertance tube is: 


dV[ 
dt 2 


A 2 


M 


( Pt~Pc ) 



The volume flow entering into volume c is described 


by: 


U - dV J V <0 d Pc 
dt yp 0 dt 


(4) 


The equation of pulse tube is: 

(5) 

YPo dt 

The length of regenerator is much shorter than the 
wavelength of working gas in the system. The momen¬ 
tum equation and continuity equation of the regenerator 
in frequency-domain are [3]: 


1 

O 

II 

(6) 

U h ~U c =F k p c+ F g U c 

(7) 

!< 

$ 

• r«o 

II 

(8) 

1-/ A 

J v g 

F k = - l -^[i+( r -i)f k ] 

(9) 

YPn, 

P _( r -!)(A-/v) 

(10) 


(1-/ v )(1-<t) 


where L n A n V r are the length, gas area of cross section, 
volume of the regenerator; cr, y are Prandtl number and 
specific heat ratio; r is the dimensionless temperature, 
which is defined as the ratio of hot end temperature T h 
and ambient end temperature T c of regenerator; f V9 f k are 
the viscous and thermal functions related to geometry 
and gas property. For ease of calculation, the randomly 
oriented stainless-steel screens are modeled by a bundle 
of parallel circular pores with the same hydraulic radius 
r h as that of the screens. In this case, the specific expres¬ 
sions of ‘f v ,fk are written as [3]: 

2J,[{i-\)2r h IS v \ 

L J 0 [(i-l)2r k /S v ][(i-l)2r k /S v ] 

2J t [(i-l)2r„ / S t ] 
k J 0 [(i-\)2rJ 8 k J_(i-\)2r h l 8 k ~] 

S v = pp/cop m 
S k =pkt(op m c p 


( 11 ) 

( 12 ) 

(13) 

(14) 


where S v , 5 k are the viscous penetration depth and thermal 


penetration depth; //, k are dynamic viscosity and thermal 
conductivity of the gas; c p is the specific heat per unit 
mass at constant pressure. 

Writing equations (6)-(7) in the time-domain, we 
have: 
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dt 


U „ 
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(15) 

(16) 


where /, X are the proportional coefficient and first-order 
differential coefficient, respectively. They can be calcu¬ 
lated by using the following complex formulation: 


z,(l + ia>A i ) = F i (i = v,k,g) (17) 


Mass conservation at the T type junction gives: 

U,=U d +U b + U R (18) 

The acoustic power dissipation outside the regenerator 
(the major part of dissipation is in the resonator) is mod¬ 
eled by an orifice and a buffer volume (which is assumed 
to be very large). The volume flow rate through the ori¬ 
fice is U b : 

U b = C 0 8p t = C'Re(/ ts ) 2 dp t (19) 
where f is the thermal function in the resonator and 

J k R 

can be calculated according to equation (12). 

The volume flow in the resonator is given by: 


U R = 


dV R V R dp 


-j- a _2_l (20) 

dt pp 0 dt 

The equation of the motion of the mass M R in the 
resonator is: 


^‘jy {sp -~ sn) <2i) 

The pressure variation in the buffer is assumed to be 
an adiabatic process and then it gives: 

Jj^dPb_ = ( 22 ) 

yp 0 dt dt 

Equations (2)-(5), (15)-(l6), (18)-(22) form an equa¬ 
tion set, from which the onset temperature and oscillation 
frequency can be known. The equation set is closed with 
eleven variables which are p h p C9 p b , U d , U C9 U h9 U t9 U b , 
Ur, V d and V R . It’s difficult to drive an analytic solution, 
so we resort to numerical calculation. 


Calculation and discussion 

In this part, some calculation results will be presented 
according to the above relations. The main dimensions of 
TASHE are listed in Table 1. The regenerator is made of 
randomly stacked screens with a porosity of 0.815 and a 
hydraulic radius of 0.25 mm. If not stated otherwise, the 
working gas is helium with mean pressure 2.0 MPa and 
the length and diameter of resonator tube are 0.35 m and 
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0.016 m respectively. 

Oscillation frequency 

For the TASHE system, oscillation frequency is an 
important parameter which is mainly determined by the 
working gas and geometrical dimensions of the compo¬ 
nents. Table 2 lists the oscillation frequencies of four 
different working gases. Calculating results of DeltaEc 
model [9] and simple lumped parameters method [2, 10] 
are also listed for comparison. It is found that the Del¬ 
taEc model predicts the frequencies with a quite high 
precision, which shows good agreement with the experi¬ 
mental results. The thermodynamic model provides a 
lower precision than the DeltaEc model but a higher ac¬ 
curacy than the simple lumped parameters method. 

The effects of resonator length on the oscillation fre¬ 
quency are shown in Fig. 2. The difference between the 
thermodynamic model and DeltaEc calculation frequency 
decreases as the length of resonator increases. When the 
length of resonator increases from 0.2 m to 0.6 m, com¬ 
pared with the experiment results, the relative error de¬ 
creases from 41.2% to 17.7%. 

Table 2 Oscillation frequencies of different working gases. 



He 

N 2 

Ar 

C0 2 

DeltaEc 

245.8 

84.0 

76.3 

67.2 

Thermodynamic 

330.0 

110.0 

103.9 

84.4 

Simple lumped 

417.3 

147.6 

135.4 

112.7 

Experimental 

245.0 

83.4 

75.5 

66.7 

Deviation (%) 

30.3 

27.5 

29.9 

26.5 


Fig. 3 gives the frequency as a function of resonator 
diameter. In small-scale TASHE, the influence of the 
resonance tube diameter on the operating frequency is 
quite sensitive. The increasing of resonator diameter will 
lead to a remarkable increase in frequency. However, the 
deviations between the DeltaEc model and thermody¬ 
namic model do not change much with resonator diame¬ 
ter. 

Onset temperature 

Onset temperature is a key parameter for evaluating 
the performance of TE. Low onset temperature has a vital 
significance on the utilization of low-grade energy [11]. 
The onset temperature of TE is influenced greatly by the 
radial dimension of regenerator, which is characterized 
by the hydraulic radius r h . A large r h means irreversible 
thermal contacts between the working gas and the stuff in 
the regenerator. If r h is big enough, the thermoacoustic 
effect will be negligibly small. 

In this case, the TE will begin self-excited oscillation 
at an extremely high temperature or even fail to self-osci- 



0.2 0.3 0,4 0 5 0.6 


Resonator length/m 

Fig. 2 Oscillation frequency as a function of resonator length 



Hate. However, the decrease of r h will lead to a small 
flow resistance of the regenerator and thus increase the 
viscous dissipation of acoustic power. Therefore, a too 
small r h will also cause a high onset temperature. There 
must be a compromise between the irreversible losses 
and the viscous dissipation of the acoustic power. Then, 
there exists an optimum value of r h for the TE to obtain 
lowest onset temperature. 

Fig. 4 shows the calculated dimensionless onset tem¬ 
perature as a function of hydraulic radius of the regen¬ 
erator. The dimensionless onset temperature is defined as 
the ratio of temperature at the hot end of regenerator to 
that at the warm end (r c =300 K). Obviously, there exists 
an optimal hydraulic radius corresponding to the lowest 
onset temperature. On the left of minimum onset tem¬ 
perature, onset temperature increases abruptly with the 
decrease of hydraulic radius, while at the right side, the 
onset temperature rises with a smaller rate. The optimal 
hydraulic radius corresponding to the lowest onset tem¬ 
perature increases with the mean pressure. Here, the 
thermal penetration depth is used to nondimensionalize 
the hydraulic radius. The dimensionless hydraulic radius 
D h is defined as the ratio of hydraulic radius to thermal 
penetration depth: 
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Fig. 4 Dimensionless onset temperature as a function of hy¬ 
draulic radius 

Equations (14) and (23) indicate that the dimen¬ 
sionless hydraulic radius varies with the mean pressure as 
well as the hydraulic radius of regenerator. Therefore, we 
can investigate the effect of dimensionless hydraulic ra¬ 
dius on the onset temperature by varying the mean pres¬ 
sure. Fig. 5(a) gives the calculated dimensionless onset 
temperatures for four different working gases as a func¬ 
tion of mean pressure. There is an optimal mean pressure 
for the TASHE to obtain the lowest onset temperature for 
each working gas. With the mean pressure increasing 
from 0.5 MPa to 3.0 MPa, the onset temperature of he¬ 
lium (He) decreases monotonically and reaches the low¬ 
est dimensionless onset temperature of 1.7 at 3.0 MPa; 
the mean pressures corresponding to the lowest tempera¬ 
ture of nitrogen (N 2 ) and argon (Ar) both are 1.2 MPa; 
the optimal mean pressure of carbon dioxide (C0 2 ) is 0.8 
MPa. When the mean pressure is less than 1.5 MPa, the 
onset temperature of C0 2 is the lowest among the four 
working gases; however, helium has the lowest tempera¬ 
ture when the mean pressure is higher than 2.5 MPa. The 
onset temperatures of Ar are always higher than N 2 from 
0.5 MPa to 3.0 MPa. 

The influence of mean pressure on the onset tempera¬ 
ture can be attributed to their effect on the dimensionless 
hydraulic radius. By calculating the dimensionless hy¬ 
draulic radius of regenerator with different mean pres¬ 
sures, we get the effects of dimensionless hydraulic ra¬ 
dius on the onset temperature as shown in Fig. 6(a). Ob¬ 
viously, the optimal dimensionless hydraulic radiuses for 
the lowest onset temperature of each working gas are 
nearly the same and the values are about 0.30-0.35. 

Pressure and Volume flow 

Proper amplitude and phase relationship of the oscil¬ 
lating pressures and volume flows have a great signifi¬ 
cance for the TASHE to obtain a good performance. Due 
to the resistance, inertial and compliance effects of the 
working gas, the pressures and volume flows in each 
component are quite different. In the thermodynamic 
model, there are three pressure variables and six volume 
flow variables as shown in Fig. 1(b). 



0.5 1.0 1.5 2.0 2.5 3.0 

Pressure (MPa) 



Fig. 5 Dimensionless onset temperature as a function of mean 
pressure: (a) in calculations; (b) in experiments. 




Dimensionless hydraulec radius 


Fig. 6 Dimensionless onset temperature as a function of di¬ 
mensionless hydraulic radius: (a) in calculations; (b) 
in experiments. 

Fig. 7 illustrates the pressure fluctuations in the time- 
domain, from which the amplitude and phase relationship 
are clearly presented. In the calculations, the pressure 
amplitude at the cold end of regenerator p c is set as 0.135 
MPa, which is a typical value from experiments. The 
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Fig. 7 Pressure fluctuations in the time-domain: (a) in calcu¬ 
lations; (b) in experiments. 

pressure amplitude at the hot end of regenerator p h is 0.12 
MPa, which has a decline comparing to p c due to the flow 
resistance effect of the regenerator. 

For the negligible inertial effect of regenerator, p c and 
p h are nearly in phase. The pressure oscillation at the 
cavity is much weaker than at the regenerator for the 
compliance effect of cavity. The amplitude of p b is 0.52 
MPa and the phase difference between p c and p b is about 
180°. 

Fig. 8 shows the volume flow fluctuations in the time- 
domain. The volume flow at the hot end of pulse tube U t 
is splited into three parts at the T type junction. The left 
one U d is fed back to the inertance tube; the right one U R 
flows into the resonator and the third part U 0 passes 
through the orifice, representing the viscous loss in the 
resonator. The relationship of the volume flows at the 
three-way tube is represented by equation (18). The phase 
difference of U d and U R is about 180°, meaning that the 



Fig. 8 Volume flow fluctuations in the time-domain. 


two flows are in opposite phase. The amplitude of U 0 is 
relatively small compared to the other three volume 
flows. 

Fig. 9 shows the phase diagram of the pressures and 
flow volumes at the two ends of the regenerator. The 
phase difference at the cold end of the regenerator is 
22.3°, with p c leading U c ; and increases to 35.9° at the 
hot end due to the compliance effect of regenerator on 
volume flow. The phase relationships calculated above 
are well matched with the DeltaEc calculation results [9]. 
The acoustic power gain, which is defined as the differ¬ 
ence between the acoustic power at the hot end of regen¬ 
erator and that at the cold end, is often used to evaluating 
the performance of the regenerator. 



Fig. 9 Phase diagram of pressure and flow volume at the two 
ends of regenerator. A quantitative scale is omitted. 

The acoustic power flowing can be written as 
E=l/2\p\\U\cos<f); therefore the acoustic power at the two 
ends of regenerator is 135.9 W and 205.1 W, respectively. 
So, the acoustic power gain is 69.2 W, which is a little 
higher than that of the DeltaEc result of 60.0 W. 

Experiments and verification 

Experiments have been performed to validate the cal¬ 
culation results. A detailed introduction to the experi¬ 
mental apparatus can be found in reference [9, 10]. The 
heating power input is fixed at 150 W. Other working 
conditions are the same as in the calculations which have 
been introduced in part 3. As shown in Fig. 1(a), pressure 
sensors labeled P c , P t and P b are located at the two ends 
of the inertance tube and the end of the buffer cavity, in 
order to compare with the three calculation pressure 
variables. The pressure sensors are of the piezoelectric 
type (CY-YD 208), supplied by Sinocera Piezotronics, 
Inc, China, with the total measuring error of ±750 Pa. A 
NiCr-NiSi thermocouples is placed in the heating block 
to measure the temperature (with accuracy ± 1°). 

The frequencies for different working gases are meas¬ 
ured as presented in Table 2. The deviations between the 
thermodynamic model calculation and experimental fre¬ 
quencies are about 30%. In essence, the thermodynamic 
model is built up based on the lumped-parameter as- 
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sumption that most components of the TASHE are 
equated to resistance, inertance or compliance. 

The treatments of these components are not accurate 
enough, especially for the small-scale TASHE in the pa¬ 
per. It has been validated indirectly that the deviation of 
frequency between the calculation and experiment de¬ 
crease as the length of resonator increases as shown in 
Fig. 2. This means the thermodynamic model has higher 
accuracy for larger TASHE system. The experimental 
values of onset temperature for different working gases 
are also presented in Fig. 5(b), 6(b). As shown in these 
figures, the calculated curves can reflect the trends of the 
experimental data. The measured onset temperature is 
higher than the computed value. The reason may lie in 
the difference between the model and the actual TASHE. 
The pressure waves are also measured as shown in Fig. 
7(b) and show good agreement with the calculated oscil¬ 
lation pressures. 

Conclusions 

The onset characteristics of a miniature TASHE were 
studied based on the thermodynamic model method, fo¬ 
cusing on the oscillation frequency and onset temperature. 
The calculation results indicate that the thermodynamic 
model can make rather remarkably accurate predictions 
about the varied trends of frequency according to work¬ 
ing gas and resonator dimensions. The onset tempera¬ 
tures of four working gases with different mean pressures 
are also calculated. Though the optimal mean pressures 
for each working gas are different, the corresponding 
dimensionless hydraulic radiuses are indeed almost the 
same, whose values lie in the range of 0.30-0.35. The 
amplitude and phase relationship of the oscillating pres¬ 
sures and volume flows are presented in the time-domain. 
The phase difference in the regenerator is close to the 
traveling-wave angle (0°), which is beneficial for the 
TASHE to obtain a good performance. Meanwhile, ex¬ 
periments have been performed to validate the calcula¬ 
tions. The frequencies and onset temperatures are meas¬ 
ured and they are in good agreement with the calcula¬ 
tions. 
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